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ABSTRACT 

Context. The contribution of major mergers to galaxy mass assembly along cosmic time is an important ingredient to the galaxy evolution scenario. 
Aims. We aim to measure the evolution of the merger rate for both luminosity/mass selected galaxy samples and investigate its dependence with 
the local environment. 

Methods. We use a sample of 10644 spectroscopically observed galaxies from the zCOSMOS redshift survey to identify pairs of galaxies destined 
to merge, using only pairs for which the velocity difference and projected separation of both components with a confirmed spectroscopic redshift 
indicate a high probability of merging. 

Results. We have identified 263 spectroscopically confirmed pairs with rj""' = 100/;"' kpc. We find that the density of mergers depends on lumi- 
nosity/mass, being higher for fainter/less massive galaxies, while the number of mergers a galaxy will experience does not depends significantly 
on its intrinsic luminosity but rather on its stellar mass. We find that the pair fraction and merger rate increase with local galaxy density, a property 
observed up to redshift z ~ 1. 

Conclusions. We find that the dependence of the merger rate on the luminosity or mass of galaxies is already present up to redshifts 
z ~ 1, and that the evolution of the volumetric merger rate of bright (massive) galaxies is relatively flat with redshift with a mean value of 
~ 3 X 10~'*(8 X 10"^ respectively) mergers M pc^^ Gyr^^ . The dependence of the merger rate with environment indicates that dense environments 
favors major merger events as can be expected from the hierarchical scenario. The environment therefore has a direct impact in shapping-up the 
mass function and its evolution and therefore plays an important role on the mass growth of galaxies along cosmic time. 
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1. Introduction 

In the last years, with large observational programs reaching ear- 
lier stages in the life of the Universe, the importance of mergers 
in the evolution of galaxies has become an observational fact. 
This is not a surprise as, in the current hierarchical structure for- 
mation paradigm, the process of coalescence of dark matter ha- 
los should strongly impact the mass assembly in galaxies living 
within these halos. 

In the hierarchical picture, the halo merging tree his- 
tory can be quantified by a halo merger rate, measuring the 
growth of mass per average mass in a representative volume 
of the Universe, and its evolution with redshift oc (1 4- z)^^^ 
( |Kitzbichler et al., 2008 ), but observational constraints come 
from the observations of galaxy-galaxy mergers which may 
not directly follow the DM merger rate evolution. In the 
last decade, a number of studies have attempted to measure 



the galaxy merger rate and its evolution using a variety of 
techniques, producing a wide range of apparently incon- 
sistent values (e.g. |Patton et al., 1997t |Patton et al., 2000| 



Le Fevre et al., 2000|Patton et al., 2002|Con selice et al., 2003 



Kartaltepe et al., 200" 



Bundy et al., 2005] 

Kampczy' ketal., 2007[ |Lotz et al., 2008t |Bundy et al., 2009| ) 
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* Based on observations obtained at the European Southern 
Observatory (ESO) Very Large Telescope (VLT), Paranal, Chile, as 
part of the Large Program 175.A-0839 (the zCOSMOS Spectroscopic 
Redshift Survey) 



but difficult to compare because of different selection functions 
and methodologies. It is only with recent large spectroscopic 
redshift surveys that robust estimate of the merger rate have 
become possible using a systematic confirmation of physically 
merging pairs, followed by a physical analysis of their proper- 
ties. In the DEEP2 survey, Lin et al. (2008) observed that the 
evolution of the merger rate has been flat with redshift. In the 
Vimos VLT Deep Survey (VVDS) sample extending to fainter 
luminosities, de Ravel et al. (2009) have demonstrated that the 
merger rate evolution strongly depends on the luminosity or 
mass of the galaxies, with brighter galaxies showing a flatter 
merger rate evolution consistent with what is observed in the 
DEEP2 sample (up to z = 1.4 at least), while fainter galaxies 
with Mb < -18 show a stronger evolution parametrized by 
(1 + 2)^^'^='"^^. In itself this result reconciles some of the 
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apparently discrepant measurements obtained previously, which 
can now be understood in the framework of galaxy samples 
with different luminosity or mass (see also Lopez-Sanjuan et al. 
(2010a)). 

In parallel, the theoretical modeling of structure forma- 
tion and large numerical simulations are produce key ad- 
vances in our understanding of the role of galaxy-galaxy merg- 
ing and associated timescales in a cosmological context (e.g. 
|Lotz et al., 2008| l. The merging timescale is an important param- 
eter when deriving the merging rate from pair count observa- 
tions, and it has been classically assumed to be roughly half a 
billion years for massive galaxies with roughly equal mass (e.g. 
IPatton et al.72000 Patton et al., 2002,Lin et al., 2004) . In a de- 
tailed analysis of the Millenium simulation Kitzbichler & White 
(2008) showed that the actual merger timescale for a pair of 
galaxies may be considerably larger for all masses than this pre- 
viously assumed value, reaching an excess of one billion years. 
Therefore, merging timescale probably remains the largest input 
of uncertainty in these kind of calculations. 

However, these simulations enable to investigate more ac- 
curately the history of galaxies with respect to their stellar 
mass accretion. There are increasingly clear evidences that the 
galaxy mass assembly process is in a very active phase be- 
tween 1 < z < 2 as observed from the evolution of the stel- 
lar mass density (e.g. Bundy et al., 2004 [Arnouts et al., 20071 
lllbert et al., 2010) , and one would naturally expect merging to 
be a major if not dominant contributor However, there is cur- 
rently a debate as to the importance of merging in the build up 
of galaxies, and in particular on how mergers impact the forma- 
tion of early-type massive galaxies at high redshift. Studies of 
the velocity fields with 3D integral field spectroscopy of bright 
z ~ 2 galaxies find a number of self supported rotating discs 
or high velocity dispersions galaxies which may be assembled 
from gas accretion via cold flows along streams of the cosmic 
web fueling the star formation rather than from merging pro- 
cesses ( |Genzel et al., 20091 |Dekel st al., 2009) 1. However, major 
merging events represent about one third of galaxies observed in 
3D spectroscopy of galaxies with z ~ 1 - 3 ([Epinat et al., 2009 
[Forster Schreiber et al., 2009| |Law et al., 2009) , indicating that 
merging is indeed very active at these epochs. 

It is clear now that mergers produce strong modifications 
of galaxies properties. The most obvious effect is a change 
in the morphology. Indeed, major mergers have been pro- 
posed to be the main mechanism that transform disk-like galax- 
ies into spheroidals. The intrinsic nature of galaxies them- 
selves could also be affected by a merging process : re- 
cent studies have shown direct evidence of merger-induced 
Active Galaxy Nuclei activity ( [Ramos- Almeida et al., 20 iT) 
[Silverman et al., 20 11) . The color of galaxies is also strongly 
affected by major merging processes. During the merger itself, 
inter-stellar clouds may collapse and many new stars are formed 
in a starbust process producing enhanced star formation rates 
(Hopkins etal., 2006; Bastian N., 2008| |de'Ravel et al., 2009l 
[Darg et al., 2009,, Bridge et al., 2010) . This enhancement of star 
formation can depletes the galaxy gas reservoir on a short 
timescale leaving the remnant "red" and "dead". We know that 
a color bimodality with a red sequence is also already present 
at z ~ 1-1.5 and that blue galaxies show only weak evolu- 
tion in number density compared to the number of red galax- 
ies that has doubled/tripled since z ~ 1 (e.g. [Ilbert et al., 2006[ 
[Arnouts et al., 2007| l. 

It therefore becomes increasingly clear that there is a chain 
of events that lower the number density of intermediate late type 
galaxies to build up more early type galaxies. The best way 



to probe the formation of galaxies that populates the red se- 
quence is to identify their progenitors. The so-called k+a galax- 
ies have been suggested to be one of them. These galaxies show 
strong Balmer absorption lines, indicative of an intense star for- 
mation epoch within the last few billion years but no-ongoing 
star formation. Multiple origins of k+a have been proposed 
(e.g. Vergani et al., 2009) but these caracteristics are those one 



could expect from a galaxy that experienced a recent wet merger 
( [Wild et al., 2009| l. Knowing the importance of wet mergers in 
the formation of k+a galaxies could therefore be a link to- 
wards the understanding of the formation of early-type galaxies. 
However, mergers also happen between galaxies within the red 



sequence. These 


are called dr}' mergers (Khochfar et al., 2003 


Bell et al., 2006b 


Di Matteo et al., 2007[ [Scarlatta et al., 2007 


Lin et al., 2008| 


Lin et al., 2010} Robaina et al., 2010i. If one 



wants to understand the role of mergers in galaxy evolu- 
tion, it is therefore crucial to understand when and where 
these different kind of mergers happen. Recent studies have 
shown that the merging activity is more active at high redshift 
( [dePropris etal., 2010| [EUche-Moral et al., 2010 ). The bulk of 
mass assembly through mergers is thought to lie between 1 < 
z<2. 



Galaxy merger rate depends on stellar mass (e.g. 
[Conselice C.J., 2006) [de Ravel et al., 20091 ) and since massive 
galaxies are more likely to live in dense regions, merger rate 
should depend on environment as well. By definition galaxy 
mergers are expected to take place in relatively high density 
regions. It has been shown, locally that , indeed, mergers 
density distribution peaks in average these intermediate/high 



density regions (e.g. Mcintosh et al., 2008 ; Darg et al., 2009 
Perez et al., 2009[ |Ellison et al., 2010 ) corresponding roughly 
to galaxy groups where the velocity dispersion is not as strong 
as in galaxy clusters. In a recent paper, Lin et al. (2010) have 
shown that this trend seems to remains at higher redshift 
(0.75 < z < 1.2). They shown a clear dependence of the 
merging rate with the local density and found that most of this 
dependence comes from the dry mergers due to the increased 
population of red galaxies in these dense environments. Our 
paper aims to give a general picture of the role of density in the 
merger rate evolution between z = 1 down to z = 0.2. 



Therefore, we report here on the use of the zCOSMOS spec- 
troscopic redshift survey of bright galaxies up to z ~ 1 to es- 
tablish a new measurement of the merger rate evolution from a 
count of real physical pairs and to test the relationship between 
the merger rate and the local environment. The layout of the pa- 
per is as follows. We present the zCOSMOS sample in Section|2] 
The methodology to identify close pairs with both members of 
the pairs having a confirmed spectroscopic redshift, and estab- 
lish the pair sample, is presented in Section |3] We detail the cor- 
rections applied in Section|4]and present the redshift evolution of 
the pair fraction in Section|5] The estimate of merger time-scales 
and the evolution of the merger rate based on the pair fraction is 
presented in Section|7] We study the influence of the luminosity 
on the merger rate in Section 16.21 Then, in Section |7] we study 
the impact of environment on the galaxy merger rate. Eventually, 
we discuss and conclude in Section|8] 



Throughout this paper we assume a flat ACDM cosmology 
with Q,„ = 0.25, Qa = 0.75 and Ho = 70 kms-^Mpc-\ All 
magnitudes are given in the AB system. 
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2. zCOSMOS Bright sample 

2.1. Spectroscopy 



zCOSMOS ( Lilly et al., 2007 1 is a large on-going spectroscopic 
redshift survey within the COSMOS dScoville et aT, 2007j field 
performed with the VIMOS spectrograph ( Le Fevre et al., 2003] l 
on the European Southern Observatory's Very Large Telescope 
(ESO-VLT). In this analysis we use the bright part of this survey, 
called the zCOS MOS-bright 10k, which covers the COSMOS 
HST/ACS field dKoekemoer et al., 2007| ). The strategy of mul- 
tiple passes of the VIMOS spectrograph will allow to reach a 
high galaxy sampling rate of 70% when complete. The 10k sam- 
ple has a mean sampling rate of 33%. Spectra of galaxies have 
been obtained using the Red Medium-Resolution grism (MR) 
with a spectral resolution R=600 covering the wavelength range 
5550 to 9650 A. This is a pure magnitude selected sample with 
Iab ^ 22.5. A total of 10644 galaxies have been observed with 
VIMOS in multi-slit mode, and the data have been processed us- 
ing the VIPGI data proccesing pipeline ( [Scodeggio et al., 2005[ ). 
A spectroscopic flag has been assigned to each galaxy provid- 
ing an estimate of the robustness of the redshift measurement 
( [Lilly et al., 20 07). If a redshift has been measured, the corre- 
sponding spectroscopic flag value can be 1, 2, 3, 4 or 9. Flag 1 
meaning that the redshift is ~ 50% secure and flag 4 that the 
redshift is ~ 99% secure. Flag 9 means that the redshift mea- 
surement rely on one single narrow emission line (Oil or Ha 
mainly). An information about the consistency between photo- 
metric and spectroscopic redshift has also been included as a 
decimal in the spectroscopic flag. In this study we select the 
highest reliable redshifts, i.e. with confidence class 4.5, 4.4, 3.5, 
3.4, 9.5, 9.3, 2.5 for main targets and for serendipitous galaxies 
within the same magnitude limit and corresponding confidence 
classes. This provides a sample of ^ 6800 galaxies in the red- 
shift range 0.2 < z < 1 with a confirmation rate of 99.5%, with 
a peak in the redshift distribution at z ~ 0.6 and a velocity accu- 
racy estimated to be ~ 100 kms^K Spectroscopic redshifts and 
photometry (see Section 12.2b have been used to compute abso- 
lute magnitude using the code ALF (see Ilbert et al. (2005) for 
further details). 



2.2. Photometry 

The input target catalogue for zCOSMOS-bright is based pri- 
marily on the F814W magnitudes derived from the HST-ACS 
images (Koekemo er et al., 2 007 ). For the small areas where the 
HST data is absent or compromised by diffraction spikes etc, 
these are supplemented by photometry from a high spatial reso- 
lution i-band CFHT image which has diffraction spikes at dif- 
ferent position angles (See Lilly et al. (2007) for further de- 
tails). The statistical reliability of the spectroscopic redshifts 
in the various Confidence Classes is assessed by the compari- 
son of repeat observations of a sample of 600 galaxies and by 
the comparison to photometric redshifts computed from codes 
taking into account stars, galaxies or AGN (Ilbert et al ., 2008[ 
ISalvato etaL, 2009; Feldma nn et al., 2008| l. The redshiffdistri- 
bution N(z) of both photometric and spectroscopic samples high- 
light two structures at z ~ 0.35 and z ~ 0.7. Moreover, in order 
to select the same regions probed with the spectroscopic sam- 
ple, we have extracted from the original parent catalogue all the 
galaxies included in the combined footprint of the VIMOS quad- 
rants. 



3. Spectroscopic pairs in tlie zCOSMOS-bright 
sample 

We identify close galaxy pairs from our spectroscopic sample by 
measuring three key quantities for each galaxy couple: the pro- 
jected physical separation r^, the velocity difference along the 
line-of-sight Av, and the absolute magnitude difference in the B- 
band AMb. Thus, close pairs are defined to have r,, < r™"' and 
Av < Av'™'. Thanks to spectroscopically confirmed redshifts, 
both projected and velocity separations between two galaxies 
are accurate enough for our study. The last criterion AM™ * , the 
maximal absolute magnitude difference, enable us to focus on 
major mergers. These particular type of mergers are those in- 
volving galaxies of comparable luminosities or stellar masses. 
As we aim to compare the merger rate of galaxies of compa- 
rable luminosities throughout the redshift range of our survey, 
we use an absolute magnitude cut of Mg = -21 at z = 1. This 
ensures that we study a complete sample at this redshift, and 
we apply a passive luminosity evolution following the evolution 
of as meas ured from the luminosity function of our sample 
( |Zucca et al., 2 009). This translates into a luminosity threshold 
evolving as M^'"(z) = Mb(z = 0) - Q x z with g = 1-36 (see 
FigD. 

Two different approaches, leading to different definitions of 
the pair fraction, have been followed. First, one can define a pair 
as a system where both members are brighter than M^"'. This ap- 
proach reduces greatly the statistics because the dynamic range 
decreases with redshift. However, this prevents from applying 
any luminosity dependent correction since both galaxies lies in 
our volume limited sample. This approach can there be seen as 
a lower limit to the pair fraction, specially at z ~ 1 where the 
reduced dynamic range reduces the observable pair statistics. 

The other way to define a pair is to apply the absolute mag- 
nitude threshold only for one galaxy, and search for companions 
satisfying the magnitude difference AMg, meaning that some of 
the companions will have an absolute magnitude fainter than the 
luminosity threshold. Therefore, this technique requires to cor- 
rect for companions we can't see in our volume limited sample at 
the highest redshift probed here. This correction relies on the lu- 
minosity function and increases with redshift since one reaches 
the limit of our volume-limited sample (see section 4.3 for more 
details). 

These two approaches answer different questions. If one 
define pairs to have both their members above the luminosity 
threshold, the pair fraction will tell us how many galaxies with 
both Mb < M^'" are in interaction. If one define pairs to have 
only one galaxy with Mb < Mg", the pair fraction will reflects 
the number of galaxies with Mb < M'™ that have a companion 
within the 1.5 magnitude range defined for the study of major 
mergers (in the B band). In the following we will present results 
using the second technique (where only one member of the pair 
has to be brighter than the luminosity threshold), however, a de- 
tailed description using the first technique will be presented in 
Kampczyk et al. (in prep). 

Looking for major mergers, i.e. equivalent luminosity /mass 
pairs with a high merging probability, imposes to use a criteria 
based on luminosity or mass difference between the two galax- 
ies. In the case of our luminosity selected sample, we use a 
criterium based on the absolute magnitude difference in the B- 
band : AM™-* - 1.5. This translates, for mass selected samples, 
in a stellar mass difference less than 1:4. Futhermore, we im- 
pose to our pairs to have a velocity difference along the line of 
sight of Av < 500 km/s. Table [1] gives the number of pairs in 
each redshift bin considered for maximal projected separations 
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of r"^"' 



20, 30, 50 and 100 h-^kpc. In the following, we will 
ignore the redshift bin < z < 0.2 because of the small vol- 
ume sampled resulting in a small number of pairs even for high 
separations. In total, 39 pairs were found across the entire field 
for r™°^ = 20 h'^kpc and 0.2 < z < 1.0. We list the properties 
of these identified close pairs with r"^-^ = 20 h^^kpc in Tabled 
This number increases to 263 pairs when we allow the separation 
to reach r™" = IQQh-^kpc. 



Table 1. Number of pairs with AM™" = 1.5 and Mb < -19.64- 
Qiz). 
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Fig. 1. The absolute magnitude in B band is plotted versus red- 
shift for all galaxies of the iOA:-sample. The solid blue line is the 
absolute magnitude cut we applied to compare galaxies of simi- 
lar luminosities throughout the redshift range of our survey. We 
assumed a passive evolution with Mb < -19.64 - Q{z)- Dotted 
lines show +1.5 magnitudes around this cut. 



4. Accounting for selection effects 

The pair fraction is basically the number of pairs over the num- 
ber of galaxies satisfying the same selection criteria. However, 
because of the selection function, estimates of pair fractions 
from spectroscopic surveys require a serie of corrections which 
are described in this section. 



4.1. Angular completeness 

We have first to investigate the angular completeness of the sur- 
vey. With a multi slits instrument like VIMOS, it is not possi- 
ble to observe two objects too close on the plane of the sky be- 
cause of the slit overlap this would produce. Furthermore, as it 
is necessary to include in each slit the sky background on each 
side of a selected target to enable for accurate sky subtraction, 
a minimum separation of about 3 arcseconds between targets is 
imposed on a single mask observation. After several VIMOS ob- 
servations of the same field using different slit masks close pairs 
are observed in a random fashion, partially erasing this minimum 
separation cut, and the fundamental limitation on pair separation 
is then set by the image quality of the images used to produce 
the photometric catalog, or about 1 arcsec. We therefore under- 
select galaxy pairs at small angular separations less than a few 
arcseconds. To estimate this bias, we compare the number of 
photometric pairs in the parent photometric catalogue with the 
number of pairs we find in the spectroscopic 70A:-sample, and we 
use the ratio between the number of spectroscopic pairs N,^ and 
the number of photometric pairs Nppas a function of the angular 
pair separation 6 as shown in Figurel2]to correct for this bias. The 
evident lack of spectroscopic pairs at small separations (< 1 ") 
corresponds to the image quality of the photometric images. For 
1 < 6 < 10 ", objects with the right magnitude can fall in the slit 
by chance, and this leads to an increase of this ratio. For 6 > 15" 
the ratio is pretty constant with a = if^iO > 15 ") = 0.0876. 

Therefore, we assign to each pair k a weight a>^ which depends 

u 

on the pair separation : 

k O 
ofk = — . 



4.2. Spectroscopic completeness 

The zCOSMOS-bright survey is not spectroscopically complete 
since about one third of galaxies have been targetted by the 
VIMOS spectrograph and measuring accurate redshifts has not 
been possible for all galaxies. This incompleteness, translated 
into a Target Sampling Rate (TSR) and a Spectroscopic Success 
Rate (SSR), must be taken into account when computing pair 
statistics and we use the same method as described in de Ravel 
et al. (2009) to correct for them. To account for the TSR, we 
apply a weight 

spectra 

_ 

lotsr - 



,jphoto 

where j^f^P^'-^'''^ jg (jjg (qj-^j number of galaxies observed in the 

spectroscopic survey and Nf'"'" is the number of galaxies in the 
parent photometric catalogue (see Lilly et al (2007) for further 
details). To correct for the SSR, we apply a weight 

spectra, flag 
^ fg_ 

where ]\f^P^''^''^^ f^'^g the number of galaxies in the spec- 
troscopic survey that have a reliable spectroscopic flag (see 
Section 12.11 ). The combined completeness weight for each 
galaxy / is then written as : 



^comp ~ 



-1 



-1 
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Table 2. List of the 39 spectroscopic pairs with rj'"" = 20h ^kpc. Ay""" - 500 km/s and AM™"" = 1.5 mag selected in the bright 
Msiz = 0) < -19.64 sample. The total number of pairs increases to 263 when r"^^' - IQQh^^kpc. 
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It depends on redshift and apparent magnitude (see 
|Zucca et al., 2009| for further details). Furthermore, we as- 
sign each pair k with 



4.3. Luminosity incompleteness 



k _ i J 

^p, comp — ^comp ^ ^comp 



where (jcomp ™^ ^comp '^^^ spectroscopic completeness 
weights of the two galaxies in pair. 

A few hundred X-ray identifications and a handful of other 
targets were included in the mask designs as compulsory tar- 
gets (see Lilly et al., 2007 for further details). Since these have 
a much greater chance of being included in the spectroscopic 
masks we include a de-weighting scheme to correct for their 
higher chance of having been observed spectroscopically. The 
mean sampling rate of these objects being 2 times that of the 
other random targets, we assigned manually an additionnal spe- 
cific weight of 0.5 to all these compulsory objects. 



Figure [H shows the distribution of absolute magnitudes in the B- 
band as a function of redshift. Here, we are looking for pairs for 
which at least one of the member is brighter than the absolute 
magnitude cut defined (see Section[3l) : Mb < Mb(z -0)-Qxz. 
At the high redshift end of our survey, we are not able to iden- 
tify pairs when the fainter companion has an apparent magnitude 
fainter than our Iab = 22.5 limit, therefore artificially lower- 
ing the number of pairs. To take this into account, we follow de 
Ravel et al. (2009) by computing for each gal axy / a weight a>'„^g 
using the luminosity function of this survey ( |Zucca et al., 2009| 
to estimate the number of galaxies thus missed. For each galaxy 
/, we derive M^^^ = + AM™' that a companion is allowed 
to have to satisfy the search criteria. Using the survey limit 
Iab = 22.5, we derive for each galaxy / at redshift z, the maxi- 
mum absolute magnitude reached by the survey at that redshift 
M^^^. We then assign a weight for each galaxy 
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Fig. 2. Completeness as a function of angular separation 0. This 
completeness is defined as the number of spectroscopic pairs A^,.. 
over the number of photometric pairs Npp. The dotted line rep- 
resent the average value at large scales with 9 > 15 ". The bump 
for separations 1 < 9 < 10 " is due to objects falling into the slit 
by chance. 



if M' ,„ < M' 



J— oo 



0(M)£/M 



J — oo 



if M' ,„ > M' 



(S>(M)dM 



where 0(M) is the global luminosity function. 
The median weight < (jj„,ag > increases as a function of redshift 
: < cJmag > = 1, 1.21, 1.58 and 2.24 within 0.2 < z < 0.4, 
0.4 < z < 0.6, 0.6 < z < 0.8 and 0.8 < z < 1.0 respec- 
tively. At this point it should be noticed that this correction is 
only based on the luminosity function with the caveat that this 
does not include any information about clustering properties of 
objects at small scales. However, Patton & Attfield (2008, Fig 
7 in their paper) have shown that, at least in the local Universe, 
at a given redshift, the co-moving number density of close com- 
panions does not depends strongly on absolute magnitude below 
over two orders of magnitude below M» . Therefore, our correc- 
tion based on the ratio between number densities of population 
with at most 1.5 magnitude difference should not introduce a 
significant bias. 

We therefore assign to each pair k of two galaxies / and j the 

■ u . k _ i j 
weignt . (i)p^ i^^^g - (^fiiag x (^jnag 



4.4. Combining weiglits 



The corrected number of pairs N™ 
fore computed as : 



and galaxies A^^"'^'^ is there- 



^ p,obs 

k = 1 



k .,..k k 

mag ^ "^g 



and the corrected number of underlying galaxies is 



N. 



g.obs 



1 



I y I 
^comp X ^mag- 



The parent galaxies from which the pairs have been identi- 
fied are those galaxies matching the absolute magnitude criterion 
within the overall iOA:-catalog. Using these corrected numbers, 
we can estimate the pair fraction fp{z) in each redshift bin as : 



fpiz) = 



Nr. 



N<^"''iz) 



Assuming that merging systems are dominated by couples 
rather than higher multiplets, which is the case in our sample, the 
pair fraction as given here represents the fraction of pairs with 
respect to the underlying galaxies. Therefore, the probability that 
an individual galaxy with Mb < M'™ has a close companion 
with AMb < AM" 



f™'" is twice this fraction. 



5. Pair fraction evolution... 

5.1. ... in luminosity selected samples 

The weighting scheme described in Section 14.41 allows to esti- 
mate the evolution of the pair fraction fp{z) with redshift. We 
perform the calculations for different sets of projected separa- 
tions (r^"*^ - 20, 30, 50 and 100/j"^ kpc) for galaxy pairs with 
^yinax _ containing at least one galaxy brighter 

than Mf'. 

To investigate the evolution of the pair fraction with red- 
shift we use the standard parametrization fpiz) - fpiz = 
0) X (1 -I- z)'". A least-squares fit, weighted by the error bars, 
to our data gives a slope m decreasing from 2.98 + 1.62 with 
r^"-* = 20 h-^kpc to 0.81 + 0.59 with r^"" = 100 h-^kpc. Pair 
fractions and corresponding best fits parameters for each sepa- 
rations are showed in Table |3] We show in Figure |3] the results 
of best fits for these different separations in the zCOSMOS sur- 
vey as well as those obtained in |de Ravel et al., 2009| within the 
VVDS-Deep sample. T he zCOSMOS-70A; and the VVDS-Deep 
( |Le Fevre et al., 2005a| l are very complementary since they have 
the same pure I-band selection function but at different depths 
(/ab < 22.5 and < 24 respectively). The VVDS-Deep sam- 
ple absolute magnitude completeness threshold is Mb < -18, 
one magnitude and half fainter than the zCOSMOS-70^. Using 
this two surveys offers the unique opportunity to extend consis- 
tently the luminosity/mass ranges of our study. At z ~ 0.9, we 
find a pair fraction for the closest pairs ir'i;"' = 20 h-^ kpc) 
of 5.00 + 1.70% for zCOSMOS galaxies whereas galaxies one 
magnitude and a half fainter (from the VVDS-Deep) show a pair 
fraction of 10.86 + 3.20% at the same redshift and for the same 
separation. Moreover, the general trend shows that the pair frac- 
tion evolution of bright galaxies (m - 2.98 + 1.62) is flatter than 
those of faint galaxies (m - 4.73 + 2.01) at any separation cri- 
terion up to z ~ 1. This indicates that faint galaxies are more 
likely to be in pairs at z ~ 1 than bright galaxies. This is not 
true anymore at z ~ 0.5 where pair fractions of both popula- 
tion are comparable. Kampczyk et al. (in prep) show that the 
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fractions of spheroid-spheroid pairs, comparable to the so-called 
t/ry-mergers, in the zCosmos-70A: sample is ~ 20% at z ~ 0.9 
whereas fraction of early-early type pairs are found to be ~ 3% 
in the VVDS-Deep. This comparaison supports the fact that the 
redshift evolution of the pair fraction up to z = 1 comes mainly 
from late type/disc galaxies since their number densities increase 
at fainter magnitudes. 

5.2. ... in mass selected samples 

In this section, we study the influence of the stel- 
lar mass on the galaxy pair fraction/merger rate. Since 
many studies have shown a clear evidence of the in- 
fluence of luminosity on the merger rate (for instance 
Patton, D. R. & Atfield, J. E , 2008] I de Ravel et al., 2009 
Lopez-Sanjuan et al., 2009a] [Lopez-Sanjuan et al., 2009b 
Bridge et al., 2010[ l, we want to understand how this depen- 
dence translates in terms of stellar masses. Stellar masses used 
for this study have been derived using the Bruzual & Chariot 
(2003) population synthesis models with an exponentially 
declining star formation history (SFH) and initial mass function 
(IMF) of Chabrier et al (2003). More details about the method 
used to derive stellar masses in the zCOSMOS-iOfc sample 
can be found in Bolzonella et al. (2009). We are then able to 
assign a stellar mass to both members of our pairs. In order to 
select major mergers, we constrain the mass ratio between the 
two galaxies to be less than 1 :4 (corresponding roughly to our 
absolute magnitude difference in the B-band of 1.5). We then 
derived for each redshift bin, the limit in stellar mass, M;,„,, 
above wich the sample is at least 95% complete for the global 
population. Depending on the lower limit M„„„ of stellar mass 
studied {logiMmmlMo) - 10, 10.5 or 1 1), we derive a correction 
based on M/,„, wich takes into account the dependence of the 
minimum detected stellar mass with redshift and mass-to-light 
ratio ( |Pozzetti eTa l., 2010). Following a similar method as 
described in section |4] we derived the pair fractions for different 
mass selected samples. In order to study the influence of stellar 
mass on the pair fraction, we had to split our sample therefore 
diminishing the average number of pairs available for each 
case. We therefore show here results for pair fractions estimated 
from pairs with r™"-^ - 100 kpc to keep enough statistics. 
Doing so, 205, 149 and 54 pairs match the stellar mass criteria 
logiM/Mo) > 10, 10.5 and 11 respectively. Figure |4] shows 
these results for different mass selected samples and compare 
them with the similar study done in the VVDS-Deep survey 
(left panel). Pair fractions values and best fits for evolution are 
shown in Table |4] 

Best fits using the standard parametrization are shown for 
the new results presented here and the VVDS-Deep. One can 
see a fairly good agreement between the two studies on the over- 
lapping selections (logiM/Mo) > 10, 10.5) except maybe for 
the highest redshift bin of logiM/Mo) > 10.5 selected galaxies. 
The zCOSMOS data show a pair fraction of / = 0.55 + 0.14 
at this redshift, higher than the one found in the VVDS-Deep 
(/ - 0.25 + 0.08) implying a steeper evolution with redshift, 
but still statistically consistent. Comparing now the two ex- 
treme selected sample, one can see a clear difference of trend 
for low-mass galaxies (logiMIMo) > 9.5) and massive galax- 
ies {logiMIMo) > 11). The pair fraction for low mass galax- 
ies raises from 0.12 + 0.04 to 0.27 + 0.05 between z ~ 0.5 and 
Z ~ 0.9 implying an evolution with m - 3.13 + 1.54 while mas- 
sive galaxies pair fraction changes from 0.28 + 0. 12 at z ~ 0.3 to 
0.20+0.07 at z ~ 0.9 implying an evolution with m = -0.14+1.2, 
consistent with a constant evolution. Therefore we conclude that 



massive galaxies are less likely to be in pairs at z ~ 1 and this 
become less and less true as cosmic time passes and the fraction 
of elliptical galaxies raises. 

Pozzetti et al. (2009) inferred from the evolution of the 
galaxy stellar mass function upper limits for the contribution of 
mergers to the mass accretion process. They found that the stel- 
lar mass accretion in intermediate mass range is mainly domi- 
nated by accretion due to star formation history. However, af- 
ter accounting for SFH mass growth they found differences of 
20 to 40% between the observed stellar mass function and the 
evolved stellar mass function. Different processes are needed to 
explain the stellar mass growth of galaxies. Among them, merg- 
ers could be a key ingredient. In that sense, these values can be 
seen as upper limits to the contribution of mergers to mass ac- 
cretion. In particular, they found a limited evolution of the stellar 
mass function above logiM/Mo) =11 which suggests that these 
galaxies formed and assembled their mass before z = 1. This 
trend is consistent with our findings that most of the merging 
activity below z = 1 appears in low mass regimes. 

6. Galaxy merger rate and stellar mass densities 
involved 

6.1. From pair fraction to merger rate 

The pair fraction gives an estimate of the number of mergers 
with respect to the underlying spectroscopic sample. Here, we 
aim to estimate the merger rate as the number of mergers per 
unit of time and comoving volume. To this effect we need to 
estimate mainly two main parameters : n(z), the comoving num- 
ber density of galaxies in each redshift bin, and r,„g ,t(r^,z*), 
the merging timescale of two galaxies with a projected sepa- 
ration of at the mean redshift z*. We take results from the 
Millennium simulations ( [Kitzbichler et al., 2008) to estimate the 
merging time-scale T,„g^ tii^p^z^), written as follows: 

T;nJ\ = ro(r*)-'/2 + X z* + /2(r*) x {logUl - 10), 

where is the stellar mass of the primary galaxy. We com- 
puted To, fx and /2 for several ranges of projected separations in 
the case of Av""" = 500 kmj s, as shown in Figure |5] (left panel). 
One can therefore assign to each pair a prediction for its merg- 
ing timescale as shown in Figure |5] (right panel). At the largest 
separations studied here (r^""^ =100 h^^kpc), the mean merging 
timescale at all redshift is ~ 3 - 3 .5 Gyr while for the closest pairs 
this timescale is ~ 1 Gyr. Estimating the merging time-scales is 
certainly the largest source of uncertainty when trying to trans- 
late pair fractions to merger rates. Depending on the method one 
use, estimates can vary by one order of magnitude (See Section 
6.3). 

Assuming that the merging time-scale is known, the merger 
rate is then expressed as : 

Nmgiz) = C,„g X fi,iz) X n(z) X r„^(z) , (1) 

where T„g(z) = 1/A^;, x 2t T„,g^ k(r^„z'') is the averaged merging 
timescale in the redshift bin considered and Cmg stands for the 
merging probability between the two galaxies. This probability 
takes into account the fact that even when spectroscopic redshift 
of both members of the pair is known, projection effects can still 
play a role. C,„g has been estimated by Patton & Attfield (2008) 
as a function of luminosity using a combination of the SDSS 
and simulated mocks extracted from the Millennium simulation. 
They found that in the local Universe, the projection effects con- 
tamination tend to increase towards fainter absolute magnitude 
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Mb < -18 - Q{z) (VVDS) Mb < -19.5 - Q{z) (zCOSMOS) 




0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 



Fig. 3. The evolution of the pair fraction with redshift for different sets of projected separations within the VVDS survey (left 
panel, de Ravel et al., 2009) and zCOSMOS (right panel, this work). Dotted lines are best fits for each separation with the standard 
parametrization fp{z) oc (1 + z)"\ 



Table 3. Pair fractions for different sets of separations and redshift with Mb < -19.64 - Q(z). 



z bin 


20h-^kpc 30h-'kpc 50h-^kpc lOOh-'kpc 


0.2 - 0.4 
0.4 - 0.6 
0.6 - 0.8 
0.8- 1.0 


0.025 ±0.010 0.047 ±0.015 0.1 12 ±0.024 0.267 ± 0.036 
0.020 ±0.009 0.062 ±0.016 0.105 ±0.022 0.297 ± 0.037 
0.057 ±0.013 0.072 ±0.014 0.140 ±0.020 0.383 ± 0.035 
0.050 ±0.017 0.093 ±0.022 0.1 19 ±0.026 0.322 ± 0.046 


Best fits for standard parametrization fp(z) = /p(0) x (1 + zY" 


fp(z = 0) 
m 


0.0088 ± 0.0072 0.030 ± 0.002 0.096 ± 0.003 0.22 ± 0.06 
2.98 ± 1.62 1.78 ±0.53 0.49 ±0.58 0.81 ±0.59 



being ~ 50% for bright galaxies and ~ 70 - 80% for galaxies 
below M,, leading to a mean value of ~ 60%. Lin et al. (2008) 
have put this constant to C,„g = 0.6 and in order to be consistent 
for comparaisons, we decided to use the same value. One has to 
keep in mind that in principle this probability should be a func- 
tion of projected separation dr, luminosity and also environment. 
Kampczyk et al. (in prep) have studied in details the impact of 
environment on the projection effects that affect pair studies. At 
the end of this section, we will estimate the range of probability 



distributions if one include a projected separation dependence in 

Here, we use C„,g - 0.6 which is the value used in de Ravel 
et al. (2009), and timescales from Kitzbichler & White (2008). 
In our analysis, the ratio C„,g/T„,g depends on the separations 
of merger candidates (both projected and radial) and on the red- 
shift. A recent study from Lin et al. (2010) have shown that this 
ratio could also evolved with the local density. However, they 
found only a weak dependancy of timescales with environment. 
On the other hand, they found that C„,g could evolve from ~ 0.8 
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Fig. 4. The evolution of the pair fraction with redshift for different mass selected samples within the VVDS survey (left, de Ravel 
et al., 2009) and zCOSMOS (right, this work). Dotted lines are best fits for each mass selection with the standard parametrization 

/p(z)<x(l+z)'". 



Table 4. Pair fractions for different mass selected samples selected with r"'" - 100 h ' kpc. 



z bin 


log(M / Ma) > 10 logiM/Mg) > 10.5 log(M/Ma) > 11 


0.2 - 0.4 
0.4 - 0.6 
0.6 - 0.8 
0.8- 1.0 


0.17 ±0.03 0.21 ±0.05 0.28 ±0.12 
0.24 ±0.04 0.22 ±0.05 0.18 ±0.08 
0.22 ± 0.03 0.30 ± 0.05 0.30 ± 0.08 
0.37 ±0.08 0.55 ±0.14 0.20 ±0.07 


Best fits for standard parametrization fp(z) = fp(0) x (1 + zY" 


fp(z = 0) 
m 


0.11 ±0.04 0.09 ±0.04 0.25 ±0.15 
1.53 ±0.74 2.29 ±0.98 -0.14 ±1.20 



in low density regions down to ~ 0.4 in the higher densities 
probed by their analysis. They found that both C,„g and T,„g are 
only marginally dependant of redshift (up to z = 1 .2). They used 
a projected separation of r™' = 50 h^^kpc to select their merger 
candidates and found a mean merging timescale of ~ IGyr. In 
our analysis, we find that for r™" = 50 h^^kpc, the typical merg- 
ing timescale is ~ l.SGyr which leads to a C„,g/T,„g ratio of 
~ 0.33 (~ 0.38 up to z=l), about half the ratio found by Lin 
et al. (Cmg/Tmg - 0.7). This comparaison show that depending 
on the way one treat the merging timescales/probabilities, there 



is a factor two uncertainty in the pair fraction / merger fraction 
translation. However, the values used in this studies are consis- 
tent with the work done on previous studies on the VVDS and 
data are therefore directly comparable. 

The previous definition of the merger rate gives then the number 
of merger events by unit of comoving volume and time, hereafter 
we will call it the volumetric merger rate. This number gives us 
an information about the global contribution of mergers to the 
mass assembly. We can look at the merger rate from a different 
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angle : the normalized merger rate. This quantity is the number 
of merger events by unit of time per galaxy. This other approach 
of the merger rate, normalized to the overall population fitting 
the selection criterium, enable us to study a potential differen- 
tial evolution for galaxies of different mass/luminosity. In the 
following, we will then give the merger rates values expressed 
from these two point of view, the volumetric and the normalized 
merger rates. 

Following de Ravel et al. (2009), we have set C„,g constant at 
0.6. Obviously, the pair fraction depends on separation but when 
translating the pair fraction to the merger rate one should get 
an absolute value. Indeed, different rj""' values should give the 
same merger rates. In Equation (1) the timescale is independent 
of C,„g and takes into account the separation effect and since 
the number density is not linked to the separations, the fluctu- 
ations we observe in the merger rate should mainly come from 
the fact that C,,,^ is put as a constant. Adding a rp dependence 
on C,ng, one should get the same merger rates for the differ- 
ent separations. A good approximation is to assume the value 
of Cingifp - 20/?"' kpc) since most timescales assumptions in 
simulations have been computed with this separation. It is then 
required that N,„g, 20 = N,ng. a for a going from 30 to lOO/i"' kpc. 
Therefore, setting C,„g{rp - 20/z"' kpc) = Cmg, 20, following 
Equation (1) the probability for a pair with a given separation 
rp - a to merge can be written as : 



f 20 X T 

mg, a ~ ^mg, 20 ^ 

^" ^ ' mg. 20 



' f-xT, 



where /" is the pair fraction with the given separation rp = a 
leading to a timescale of T^g^ „ . We show in Figure |6]the evolu- 
tion of Cmg, a with C,„g, 20- 

These estimates give us a range of merging probabilities as 
a function of rp and z. For instance, if one put C,„g^ 20 = 0.6, 
the merging probability for a pair with similar stellar masses is 
0.30 -0.70 at = 30 h'^ kpc, 0.25 -0.55 at r,, = 50 /i^'/tpc and 
0.15-0.4 at r^, = 100/i 'A;/7c depending on redshift. Fluctuations 
seen in merger rates could therefore be explained by introducing 
a dependence in Cmg with separation and redshift. The a poste- 
riori values we found for C,„g{rp, z) are reasonnable and follow 
a consitent trend. This ensure us that, at least qualitatively, the 
numbers we derived using Millennium simulations are consis- 
tent. 



6.2. Influence of luminosity on the merger rate 

We investigate here the dependency of the evolution of the 
merger rate on luminosity. We give values of volumetric merger 
rates for the luminosity selected sample in Table |5] In order to 
quantify the evolution with redshift of this merger rate, we use 
the standard parametrization N,„g(z) = Nmg(0) x (1 -n z)™"*. We 
performed a least-squares fit to our data and best fit parameters 
are also given in Table |5] Depending on separation, we found 



that the evolution index m^g evolves from (1 + z) 



2.79±2.14 



with 



^ax ^ 20 h-^kpc to (1 + 2)0-27±0.92 ^jjjj ^ax ^ ^qq j^-lj^p^ 

Since merger rate should not depend on the projected separation, 
we can combine the different best fit values found for the differ- 
ent separations and estimate that the mean merger rate evolves 
as (2.04 + 0.83)(l -Hz)°''^=^° '*2 (xlQ-^ mergers Mpc'^ Gyr'^). 

Normalized merger rates (i.e. the number of mergers per unit of 
time per galaxy) at different redshifts are also given in Table |6] 

We show in Figure |7] the evolution of both volumetric and 
normalized merger rates for various projected separations in the 



Table 5. Volumetric merger rates for different sets of separa- 
tions and redshift with Mb ^ 
Mpc^^ Gyr'^. 



< Mi'" in units of xlO mergers 



z bin 


20h-'kpc 30h-'kpc 50h-'kpc lOOh-^kpc 


0.2 - 0.4 
0.4 - 0.6 
0.6-0.8 
0.8-1.0 


1.87 ±0.46 2.43 ±0.46 3.74 ± 0.48 4.63 ± 0.38 
1.05 ±0.28 2.24 ±0.36 2.48 ±0.31 3.69 ± 0.28 
3.39 ±0.46 2.90 ±0.34 3.69 ± 0.32 5.37 ± 0.29 
2.73 ± 0.54 3.42 ± 0.49 2.91 ± 0.38 4.22 ± 0.36 


Best fits for N,„g{z) = N„,i,{0) x (1 + z)"""' (xlO"") 


ffll})g 


0.55 ±0.59 1.53 ±0.34 3.16 ± 1.65 3.92 ±1.75 
2.79 ±2.14 1 .23 ± 0.72 -0.03 ± 1 .07 0.27 ± 0.92 



Table 6. Normalized merger rates for different sets of 
separations and redshift with Mb < M^'" in units of 
xlQ^^mergers Gyr^^ gcil^ ■ 



z bin 



0.2 - 0.4 
0.4 - 0.6 
0.6 - 0.8 
0.8 - 1.0 



20h-'kpc 30h-'kpc 50h-'kpc imr^kpc 



1.98 ± 0.48 2.56 ± 0.48 3.95 : 

1.55 ±0.41 3.30 ±0.52 3.66: 

4.43 ± 0.61 3.79 ± 0.45 4.82 : 

3.79 ± 0.75 4.75 ± 0.69 4.04 : 



0.51 4.89 ±0.40 

0.46 5.44 ± 0.41 

0.42 7.01 ± 0.38 

0.52 5.86 ± 0.50 



zCosmos iOfc-sample for galaxies with Mb < -19.64 - Q(z). 
We compare the zCOSMOS results to measurements derived 
from the VVDS-Deep sample ( |de Ravel et al., 2009| l. Within 
the iOA:-sample (Iab ^ 22.5), the mean merger rate of bright 
galaxies evolves as (2.04 + 0.83)(1 + 2)0''v±0-82 (xiQ-^ merg- 
ers Mpc^^ Gyr^^). We estimated in the same way the mean 
merger rate derived within the VVDS-Deep (Iab ^ 24) and 
found that it evolves as (4.82 + 1.62)(1 +2)206±o.62 (xiQ-^ merg- 
ers Mpc^^ Gyr^^) for faint galaxies. Comparing these results, 
we find that the evolution of the merger rate is steeper in fainter 
samples. This confirms and extends to brighter luminosities the 
trend seen in |de Ravel et al., 2009| Integrating these volumetric 
merger rates over time between z = and z = 1 is a way to quan- 
tify the contribution of different galaxy populations to the global 
merger history in the Universe. We find that galaxies brighter 
than Mb = -18 - Q(z) contribute ~ 3.8 times more to the global 
merger history than galaxies brighter than Mb = -19.6 - Q{z) 
up to z = 1 . 

Figure |7](Wg/!f panel) shows the normalized merger rates per 
galaxy in the two samples. The first thing one can notice is that 
the number of mergers per galaxy has decreased by a factor of 
2 between z ~ 0.9 and z ~ 0.3, evolving from about ~ 0.06 to 
~ 0.03 respectively. We don't see significant effects of the lu- 
minosity on this evolution. The faint galaxies seem to show a 
higher normalized merger rate but this trend, according to the 
size of error bars (reflecting poissonnian errors), does not have 
a high statiscal significance. This means that a given galaxy, re- 
gardless of its luminosity, will experiment two times less merg- 
ing processes between z = 1 and z = 0.3. 
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Fig. 5. Left Panel : Evolution of merging timescales with redshift for different projected separations rp. Right panel : Merging 
timescales of galaxy pairs as a function of separation rp for different redshift bins. 



r™"^ = 30 h~^kpc r"^"'' = 50 h'^kpc r^"^' = 100 h-^kpc 




Fig. 6. Probability of merging as a function of Cmg, 20 for different redhift bins. Panels from left to right give the probability to merge 
for various separations (30, 50 and 100 h'^kpc) as a function of Cmg, 20 

. The black dotted line is the same in each panel and is just plotted as an eye-guide. 



6.3. Influence of stellar mass on the merger rate 

Figure |8] shows the volumetric galaxy merger rate as a function 
of redshift for the different mass-selected samples. For interme- 
diate masses {log{M jM^) > 10 and 10.5) we find good agree- 
ments with results from the VVDS. This enables us to compare 
the lowest (logiM/Me) > 9.5) and highest (logiM/Me) > 11) 
mass regimes directly. 

There appears a clear difference of volumetric merger rate 
trends between high and low mass galaxies. We find that low 
mass galaxies {log{M IMq) > 9.5) are more likely to have a ma- 
jor companion than high mass galaxies (logiM/Mo) > 11) by 



a factor of order of ~ 10 below z = 1. For the less massive 
sample in the VVDS {logiMIMo) > 9.5), we find a merger rate 
evolving as Nmg(0) x (1 + z)"'™? with m,„g - 2.38 + 1.57 and 
N„,g(z = 0) = (3.56 ± 3.17) x 10"^ mergers h^Mpc-^Gyr-\ 
while for the massive sample in zCOSMOS-iOA; sample, we 
find that the merger rate evolves with m,„g = 0.73 + 1 .80 and 
N„,g(z = 0) = (0.48 + 0.44) x 10""^ mergers h^Mpc-^Gyr-^ 

Therefore, major mergers below z = 1 are found to be more 
likely for low mass galaxies, merging history of high mass galax- 
ies being roughly low and constant for massive galaxies during 
the last eight billion years of galaxies evolution. We find an av- 
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Fig. 7. Left panel : Evolution of the volumetric galaxy merger rate with redshift for different projected separations in the zCosmos 
70^-sample (solid lines). We show in dotted lines the merger rates derived from the VVDS-Deep survey (Iab ^ 24) one magnitude 
and half deeper.Right panel : Evolution of the merger rate per galaxy in the zCosmos iOA:-sample (solid lines). Dotted lines show 
again the comparaison with the VVDS-Deep survey. 



erage value of 8 x 10"^ mergers h^Mpc^^Gyr^^ for our massive 
sample in really good agreement with the results found by Chou 
et al. (2010) on the flat evolution of massive dry mergers up to 
z = 0.7. 

Figure |8] also shows a comparison between our estimates of 
volumetric merger rates and several previous studies which use 
different methods to infer the galaxy merger rate. 

Kitzbichler & White (2008) made N-body simulations 
of DM evolution (Millennium) in combinaison with a semi- 
analytical model for galaxy formation to explore the relationship 
between the evolution of galaxy merger rate and observed close 
pairs. They find that the galaxy merger rate decreases with stel- 
lar mass. We find good agreements with their data, except maybe 
for the most massive galaxies sample, where we find a moder- 
ate ( more or less flat ) evolution with redshift while they find a 
strong decrease. 

Bundy et al. (2009) used deep Infra Red observations on 
the northern GOODS field combined with public surveys in 
GOODS-S to investigate the dependence on stellar mass of the 
merger rate up to z ~ 1 .2. At least up to z = 1, the highest red- 
shift we can reach here, merger rates are in good agreement. 

Lopez-Sanjuan et al. (2009a) measured the fraction of galax- 
ies undergoing disk-disk major mergers at 0.35 < z < 0.85 
by studying the asymmetry index A of galaxy images. They 
have corrected with care the bias due to varying spatial reso- 
lution and image depth with redshift and developed a method 
to account for incertainties in asymmetry indices and photomet- 
ric redshift. Morphological and spectroscopical studies of close 
pairs are hard to compare since, by definition, the two technics 
trace mergers at different stages. In the litterature, morphologi- 
cal studies tend to find higher pair fractions than those derived 
with spectroscopic samples. Nevertheless, we find good agree- 
ment with Lopez-Sanjuan et al. (2009a) which can be seen here 



as a lower limit to the global merger rate, since they studied disk- 
disk mergers. We also compare our results with Lopez-Sanjuan 
et al. (2009b). In this paper, they measured the merger rate of 
log{M) > 10 galaxies in GOODS-S using also assymetries. 
Again a nice agreement is found. 

In another morphological based study, Conselice et al. 
(2009) used deep HST imaging from the Extended Groth Strip 
and the COSMOS field to examine structural properties (CAS) 
of galaxy pairs and derived the merger rate evolution for galax- 
ies with logiMIMo) > 10 up to z = 1.2. Their results are found 
to be in disagreement with ours by a factor of 2-3 at all red- 
shifts considered here whereas pair fractions seem consistent. 
Lotz et al. (2008) showed that the time-scale for a 20 kpc pair 
to merge is 0.2 + 0.1 Gyr, significantly lower than what we 
used in this analysis («; \Gyr). This might be the main expla- 
nation for the apparent disprecancy we find in merger rates. 
Conselice et al. (2009) also investigated the systematic differ- 
ence between morphological-based and pair count estimates of 
merger rates. They compared pair fractions derived with pair 
count (|Kartaltepe et al., 2007J and morphological criteria and 
found that this apparent disprecancy (by a factor varying from 
1.5 to 3.5 with Kartaltepe et al. (2007) results) comes from the 
different time-scales probed by the two methods. However, al- 
though they used the same approach, our agreement is much 
better with |Lopez-Sanjuan et al., 2009b raising up an issue be- 
tween morphologically based analysis themselves. Lopez et al. 
(2009b) claim that this discrepancy with Conselice et al. (2009) 
is due to the fact that they don't take into account the obser- 
vationnal errors contained in the asymmetry parameter A (see 
Section 6 of their paper). This effect could lead to an over- 
estimation by a factor ~ 3 of the final merging rate consistent 
with what the comparaison shows. 
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Figure |9] shows the evolution of the normalized merger rate 
for different mass selected samples in the zCOSMOS and the 
VVDS surveys. As for the luminosity selected samples, one 
can see a weakening of the trend when one normalizes to the 
global population. However this effect is less strong for mass 
selected samples since massive galaxies {log{M IMq) > 11) 
show between two and three times less merging activity (per 
galaxy) than less massive galaxies (logiM/Mo) > 9.5) at z ~ 1. 
Indeed, Kitzbichler & White merging time-scales are found to 
enhance the merger rate stellar mass dependence while their de- 
pendence with luminosity is less strong (Bundy et al., 2009). 
This could explain the different trends we see here between stel- 
lar mass and luminosity selected samples. Anyway, this result 
remains statistically significant and shows that low mass galax- 
ies (logiM/Mo) > 9.5) are more likely to experience a major 
merger at high redshift. Moreover, at any redshift below z = 1 , 
these galaxies show more merging activity than the most massive 
galaxies in our sample {logiMIMo) > 1 1). 

6.4. Stellar mass density involved in merger processes 

The volumetric merger rate estimates the number of galaxy 
mergers by unit of time and volume. Assuming that one know 
the total stellar mass involved in the merging process, this can be 
translated into the fraction of the stellar mass which is involved 
in a merger process. For each pair k, containing galaxies / and j, 

we are able to define the stellar mass M^^^g^^ - Mi + Mj in- 
volved in the merger process. The stellar mass density involved 
in the merger is then written as : 



merger = Nmgiz) X r^,-„(z) X 



Np{z) 



Np(z) 
k^O 



M 



merger 



(z), 



where Tbm(z) is the elapsed time correponding to the redshift bin 
concerned, A^^ the number of pairs. We estimate the fraction fp 
of the stellar mass density involved in a merging process relative 
to the underlying population as : 



merger, , 



where P]\^(z) is the mean stellar mass density of the global pop- 
ulation for each redshift bin. 

Following Patton et al. (2000), we estimate the fraction of 
stellar mass that have undergone a major merger since z ~ 1 . We 
performed this calculations for mass selected galaxies and found 
that 12.2 + 3.9, 13.1 + 5.6 and 9.7 + 5.2% of stellar mass with 
logiM/Mo) > 10, 10.5 and 11 has been involved in a merger 
process since z = 1. This shows the higher merging activity 
trend of less massive galaxies below z = 1 . Given the uncer- 
tainties on the merging time-scales and the stellar masses them- 
selves, we definitely can't rule out a consistent merging history 
between those galaxy populations at this point. However, these 
numbers can be compared to the number of ~ 20% for galaxies 
with logiM/Mo) > 9.5 found in de Ravel et al. (2009), show- 
ing that the major merging activity for today massive galaxies 
(logiM/Mo) > 11) since z = 1 has two times less important 
than lower mass galaxies (logiM/Mo) > 9.5) suggesting an ear- 
lier epoch for massive galaxies mass assembly through merging 
process. Out of the 20 to 40% of residual evolution found by 
[Pozzetti et al., 20 1"0| after accounting for star formation history, 
between 10 and 15% are due to major mergers. In other words. 



between z ~ 1 and z ~ 0, our studies suggest that stellar mass ac- 
cretion of galaxies is due to in situ star formation activity at the 
level of ~ 70% and major mergers at the level of ~ 10 - 15%. 
This leaves 15 - 20% of evolution that could come from other 
physical processes like minor mergers or prolongated star forma- 
tion histories (second bursts). Lopez-Sanjuan et al. (2010b) have 
shown that the relative contribution of the mass growth by merg- 
ing is ~ 75% / 25% due to major / minor mergers. In our case 
this would translate into a minor merger induced mass accretion 
of ~ 2.5 - 4% since z ~ 1. 



7. Impact of environment on the global merger rate 

The galaxy overdensity field in the zCOSMOS area has been 
studied in details in Kovac et al. (2010). Using this density infor- 
mation, we are able to investigate in which type of environment 
mergers are preferentially taking place as cosmic time evolves. 
Many density tracers have been derived to estimate the density, 
and different tracers of galaxy populations have been used. Here 
we used the flux limited sample {Iab ^ 22.5) which enable us to 
use our entire sample in order to keep enough statistics. The den- 
sity computations have been performed either on a fixed scale ( 
for instance 3 Mpc/h) or on variable scales corresponding to the 
projected distance of the N''' nearest neighbour The mean sep- 
aration between a galaxy and its 5'* nearest neighbour has been 
shown to be always less than 3 Mpc/h. This approach traces the 
highest densities of the field, with the caveat that the scale varies 
with redshift (from ~ 1 Mpc/h at z ~ 0.2 to ~ 3 Mpc/h at 
z ~ 1). Since we are interested in the evolution of the mergerrate 
with redshift, we decided to use densities computed using a fixed 
aperture of 3 Mpc/h. Densities have been computed either using 
a simple count of galaxies, or using a weight on the luminosity or 
the stellar mass of galaxies. It is expected that some galaxy prop- 
erties, and particularly masses, are better tracers of the underly- 
ing matter density field, therefore we decided to use the densi- 
ties computed using a mass weighted scheme. We studied the 
influence of the weighting scheme on our results and found that 
the observed trends remain regardless of the weighting scheme 
used. Therefore we decided to use densities estimated with the 
flux limited sample, on a fixed scale of 3 Mpc/h and where the 
density is weighted by the stellar masses of galaxies. In the gen- 
eral case, the density at a given position r = r{a, 6, z) is given by 
dKovacetal., 20101 ): 



Pir) = 2 



miW{\r-ri\-R) 

4>(n) 



where is a spatial window function, m, is the weight based 
on the stellar mass of galaxies, 0(r,) a correction function for 
various observationnal issues and R is the smoothing length. In 
the following, we will refer to the density (density contrast) 6 
defined as : 

. Pir)-Pmiz) 
S(r) = , 

Pm(z) 

where 0m(z) is the mean density at a given redshift. 

Figure [To] shows the density distribution of the underlying sam- 
ple and the pair samples (selected with rj""^ = lOO/i"' kpc and 
Mb < -19.64 - Q(z)) within several redshift bins. Galaxy pairs 
are found to be systematically in denser regions at all redshifts. 
This is of course expected but to quantify this trend, we studied 
the evolution of the fraction of pairs with respect to the underly- 
ing sample as a function the density 1+6. 
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Fig. 8. Evolution of the volumetric merger rate as a function of redshift for different mass-selected samples from zCosmos (red 
filled circles) and VVDS (blue empty circles). We also plot several previous studies in black: Conselice et al.,2009 (empty circles), 
Lopez-Sanjuan et al., 2009a (empty triangle), Lopez-Sanjuan et al., 2009b (filled squares). Bell et al., 2006a (cross), Bundy et al., 
2009 (empty squares) and Kitzbichler & White, 2008 (dotted Une) 
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Fig. 9. Evolution of the merger rate per galaxy as a function of redshift for different mass-selected samples from zCosmos (red filled 
circles) and VVDS (blue empty circles). 
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Fig. 10. Over-density distribution of the underlying sample (top 
panels) and r""^ - lOOh^^kpc selected pairs (bottom panel) 
for several redshift bins. We used fixed aperture {3Mpc) mass 
weighted densities estimations. The dashed line is the median of 
the distribution in each case. For all redshift bins, the distribution 
median is slightly higher for paired galaxies. 



Figure [TT] shows the pair fraction as a function of the over- 
density 6 for 0.2 < z < 0.5 and 0.5 < z < 1 and for different 
projected separations r™"-^ - 100/;"' kpc. 
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Fig. 11. The fraction of pairs with respect to the underlying 
sample is plotted versus the over-density for pairs with rj™ ' - 
lOO/z ' kpc. We show in dotted circles the shift applied to the 
low redshift sample to correct for the evolution of the mean den- 
sity of the structures (see text). The dotted line is the best fit we 
get while fitting the high-redshift and corrected low redhift bins 
data. 



Abscisse values are defined using quartiles of densities dis- 
tributions. We find that the fraction of pairs is higher in denser 
environments for both high and low redshift regimes. In order 
to evaluate the evolution with redshift of this relation, we have 



corrected the densities from the evolution in the mean density 5m 
evolving as oc (1 -Hz)-', translating to a change in the mean density 
of A(Zog( 1 + 6m)) ~ 0.45 between the two redshift ranges consid- 
ered, as plotted in Figure [TT] Comparing the 'mean density cor- 
rected' relations at low redshift and high redshift values, we in- 
fer that the fraction of pairs depends mainly on the local density 
with little or no evolution with redshift, consistent with recent re- 
sults within the DEEP2 survey by Lin et al. (2010). A trend also 
found in samples with less statistics with r™°-^ - 30, 50 h^^ kpc. 
The dependence of the number of pairs to the local density, over 
the full redshift range considered and up to z ~ 1 is compati- 
ble with a power law evolving like (1 + (5)~'' '^, (1 + 6)~^^^^ and 
(1 -H(5)~" ""' for separations of 100, 50 and 30 h^^kpc, respectively. 
This means that for an environment 10 times denser, one finds 
3.0, 2.6 and 2.1 times more pairs for separations of r""* - 100, 
50 and 30 h^^kpc respectively. Finally, quantitatively speaking, 
the pair fraction is about 2-3 times higher in environment one 
order of magnitude denser 

Kampczyk et al. (in prep) have investigated the possible con- 
tamination due to objects that fit our pair criteria but will not 
lead to a merger process, which corresponds to the introduction 
of a environment dependence in C,„g. This effect is expected to 
be stronger in over-dense regions since the gravitationnal bound 
between objects is weaker They find that for the highest over- 
density quartile this correction can decrease the pair fraction by 
20% for pairs with dr'™^ - 100 Ir^kpc and that it is consid- 
erably smaller for lower separations of t/r""" < 50 h~^kpc (up 
to ~ 12%) and t/r'™* < 30 h'^kpc (less than 4%). For all the 
other over-density quartile, this effect is almost negligible (see 
also Lin et al. 2010 for a theoritical perspective). 

In order to estimate the contribution of different density quar- 
tiles to the global merger rate we show in Figure [T2]the merger 
rate as a function of environment for different redshifts bins us- 
ing < 100 /i"' kpc pairs. We believe this estimate is legit be- 
cause our merger rate estimation takes into account the fact that 
< 100 /i"' kpc pairs have higher merging time-scales on aver- 
age and in theory we should be able to get a merger rate re- 
gardless of the separation assumed to define pairs. As expected, 
most of the global merger rate comes from the high density re- 
gions. However, the impact of local density on the merging time- 
scale is not taken into account here but, as mentionned before, 
this effect has been estimated to be up to ~ 20% for the very 
dense regions. This will lead to a decrease of our merger rate es- 
timate from ~ 2 to ~ 1.6 x lO""* mergers h^Mpc^^Gyr^^ in the 
0.8 < log(l + d)< 1.0 regions. 

We find the evolution of the merger rate with environment to 
be weakly dependant on the redshift up to z = 1 . We performed 
a least square fit on data to parametrize this evolution. We found 
that the galaxy merger rate evolves as N,„g(6) - ax logi \+6) + b 
with a = (1.61 + 0.13) X lO^* and b = (0.41 + 0.05) x lO"". 
This means that the number density of merging processes in high 
density regions is increased by a factor ~ 8 - 10. 



8. Discussion and conclusions 

We have studied the evolution of the merger rate with redshift 
in the zCosmos-70fe bright survey for galaxies brighter than 
Mb < -19.64 - Q{z). We have identified a sample of 39 (263) 
pairs with each member of the pair having a secure spectroscopic 
redshift with separations up to (20(100) h^^kpc, 500 km i"'). 
With the avaibility of the local galaxy density field measured for 
this sample, we have been able to investigate the influence of 
environment on the pair fraction and the merger rate. The pair 
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Fig. 12. Merger rate as a function of environment for different 
redshift bins. The dotted indicate the best for all the data and the 
shaded area the 1 - cr error on the best fit. 



fraction increases from 2.5 + 1.0% at z ~ 0.3 to about 5.0 ± 1.7% 
at z ~ 0.9 for galaxies selected with Iab < 22.5 and brighter than 
Mb = -19.6 - Qiz). When comparing to galaxies also 1-band se- 
lected in the VVDS-Deep survey, but 1.5 magnitude fainter, with 
a pair fraction of about 10% at z ~ 0.9. (de Ravel et al., 2009) , 
we find that the pair fraction is dependent on luminosity, with 
bright galaxies showing a lower pair fraction, confirming the de 
Ravel et al. (2009) result. We find that the pair fraction of bright 
galaxies evolves with redshift as (1 -H z)'" with m - 2.98 + 1.62 
for separations of 20 /i"' kpc decreasing to m = 0.81 + 0.59 
for separations of 100 /i ' kpc. Using merging timescales de- 
rived from the Millennium simulation dKitzbichler et al.,2008| l, 
this translates to a mean merger rate evolving as (2.04 + 
0.83) X (1 + 2)0-97±o.82 (xio-4 mergers h^Mpc'^Gyr-^) 
for galaxies brighter than Mb < -19.64 - Qiz). The mean 
merger rate between z = and z = 1 is about 3 x 
lO"'* mergers h^Mpc-^Gyr-^). This is ~ 5 times lower 
than what has been found for galaxies 1.5 magnitude fainter 
( |de Ravel et al., 20"09| . We therefore confirm that the pair frac- 
tion and the volumetric merger rate are dependant on luminos- 
ity, with fainter galaxies involved in more mergers than brighter 
galaxies. This results is consistent with what is found in the 
analysis of the Millenium simulation (Kitzbichler et al ~2008[ 
IPatton, D. Rr& Atfiel d^ J. E., 2008| l. However, when normal- 
ized to the global underlying population, the number of merg- 
ers per galaxy of both bright and faint population are compara- 
ble while less massive galaxies still show 2-3 times more merg- 
ing activity than massive galaxies at z ~ 1 . This indicates that 
the density of mergers depends on the depth of the observa- 
tions, being higher for fainter observations, while the number 
of mergers a galaxy will experience does not depends strongly 
on its intrinsic luminosity but rather on its stellar mass. This re- 
sult indicate that the merging time-scales derived by Kitzbichler 
& White (2008) are more sensitive to the stellar mass. This is 
confirmed when comparing these time-scales to those estimated 



from Patton & Atfield (2008) where the stellar mass dependence 
of time-scales is shown to be more moderate. 

We found a clear difference of merging histories between 
low and high mass galaxies, volumetric merging rates evolving 
as A^;„^(0) X (1 + z)'"""^ with m,ng = 2.38 + 1.57 and N„g{z = 
0) = (3.56 + 3.17) X 10""* mergers h^Mpc-^Gyr-^ for galax- 
ies with (logiM/Mo) > 9.5), while for the massive sample in 
zCOSMOS-iOfe sample, we find that the merger rate evolves 



0.73 + 1.80 and NUz = 0) = (0.48 + 0.44) x 10" 



with m, ^ 

mergers h^Mpc-^Gyr-^ . This leads to a fraction of ~ 20% 
(~ 9.7%) of the today stellar mass with logiM/Mo) > 9.5 (re- 
spectively > 1 1) that have been accreted through a major merger 
since z = 1. 

We have measured the pair fraction as a function of the local 
environment and we find that galaxy pairs are more likely to be 
found in high density regions with a probability increasing by a 
factor of 2-3 when increasing the local density by a factor of 10. 
The pair fraction is globaly higher at higher redshift at all densi- 
ties, following the evolution of the mean density of the Universe 
with cosmic time. As expected, the pair fraction increases faster 
in high density regions when the pair separation is increased. 
We have derived the evolution of the merger rate as a function 
of environment and find that the merger rate is increasing by one 
order of magnitude when the local density increases by the same 
amount, parametrized with Nmg{6) -ax log(l + 6) + b with 
a = (1.61 + 0.13) X 10-^ and b = (0.41 + 0.05) x lO^^^. We 
find that the values of this parametrization are only marginally 
different as redshift increases. Recent local measurements of the 
evolution of the fraction of mergers as a function of environment 
have shown that the peak of mergers distribution resides in inter- 
mediate density regions ( |Darg et al., 2009 1. These results seem 
to be qualititavely in agreement with our conclusions on higher 
redshifts. From these measurements, we therefore conclude that 
higher environments produce higher merger rates with only a 
weak dependence with redshift up to z = 1 . This indicates that 
the environment has a clear and constant effect on the merger rate 
throughout cosmic time. Galaxies in higher density regions are 
subject to 10 times more merging events when the density is in- 
creasing by a factor of 10 whatever the epoch between z = 1 and 
z = 0. This behaviour is naively expected as in higher density re- 
gions the mean separation between galaxies decreases. The link 
between merger rate and local density is also supported by simu- 
lations ( |Perez et al., 2009| l. Galaxies in high density regions will 
therefore accrete mass significantly faster than in low density 
regions. This is expected to significantly contribute to several 
properties like the observed differences in the stellar mass func- 
tion and luminosity function of galaxies in high and low density 
environments (Bolzone lla et al., 2 010; Zu cca et al., 2 009). This 
study therefore brings new evidence that the environment plays 
a key role in shapping up galaxies along cosmic time. 
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